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Abstract 

We discuss the energy level splitting Ae due to quantum tunneling between congru- 
ent tori in phase space. In analytic cases, it is well known that Ae decays faster 
than power of h in the semi-classical limit. This is not true in non-smooth cases, 
specifically, when the tori are connected by line on which the Hamiltonian is not 
smooth. Under the assumption that the non-smoothness depends only upon the 
X- or p-coordinate, the leading term in the semi-classical expansion of Ae is de- 
rived, which shows that Ae decays as fi^~^^ when ?l — > with k being the order of 
non-smoothness . 
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1 Introduction 



This paper concerns the sphtting of classicahy degenerate energy levels. The near degeneracy 
(ND) classically corresponds to congruent tori in phase space while the quantum tunneling 
between the tori causes the splitting A well-known example occurs in the one-dimensional 
symmetric double-well potential, where the eigenenergies below the top of the barrier cluster 
into two- fold ND's with energy difference vanishes as 

Ae ~ ?i"e-^/^ (1.1) 

when h ^ 0. When turning to multi-dimensional cases, M. Wilkinson showed that Ae vanishes 
normally in the same or, in certain situation, even more singular manner than (1.1) However, 
is it always true that the energy level splitting resulting from quantum tunneling is smaller than 
any power of h in the semi-classical limit? Let us see the following example. 

Consider the system on one-dimensional circle defined by any of the four Hamiltonians, 
ifi = — + cos^ x, ii'2 = + I cosx|, = \p\ + cos'^ X , H4 = \p\ + \cosx\. (1.2) 

In classical mechanics, the above Hamiltonians determine similar phase space portraits, particu- 
larly, the motion at H ^ 1 contains two symmetric closed orbits: the two vibrational orbits with 
H < 1 are connected by the transition {x,p) (x + iT,p) and the two rotational orbits with 
H > 1 are connected by the time reversal {x,p) (x, —p). According to the Einstein-Brillouion- 
Keller (EBK) semi-classical quantization rule, this classical degeneracy implies a two-fold ND 
structure in the spectrum of H. We can verify this prediction by directly diagonalizing the 
Hamiltonians. In Fig.l, Ae is plotted as the function of the mean energy of the ND pair(e). 
As expected, Ae (open dots) is much smaller than the spacing of e (approximately the dotted 
lines). However, contrary to the exponentially decay of Ae with the increase of |e — 1| in Fig. 1(a), 
non-exponentially decay of Ae in some cases is obvious. From the four illustrations, we can see 
that the "exceptional" ND occurs when and only when the corresponding classically degenerate 
tori (closed orbits) in phase space are connected by line(s) where the Hamiltonian is not smooth. 
This fact suggests that tunneling between the degenerate tori can be greatly enhanced by the 
passage of non-smoothness. 

In fact, M.V. Berry showed this non-smoothness-enhanced quantum transition between clas- 
sically degenerate states about two decades ago[^. In studying the coefficient r for reflection 

^In this paper, the word "tunneUng" refers to the quantum transition between states that classically correspond 
to separate tori in phase space M. 
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above a barrier V{x) in the semi-classical limit, Berry proved that r ~ when V{x) has a 
discontinuous k-th derivatives, in contrast to the analytic case where r is exponentially small. 
Another interesting quantum manifestation of non-smoothness, the power-law localization of 
eigenstates was also discussed in more recent papers(e.g. [^, |^]). 

In this paper we shall investigate the energy level splitting resulting from the non-smoothness- 
enhanced tunneling. We first consider the case where ND is related to the time reversal symme- 
try. By perturbation method, a relation between Ae and the non-smoothness of the potential is 
derived. Based on a geometrical interpretation, this relation is applied to a class of non-smooth 
systems. 



2 Power-Law Energy Splitting 



In this section we study systems where ND is related to the time reversal symmetry. The 
problem is more tractable since the projection of torus onto the coordinate space contains no 
singularity (caustic). By perturbation method, we obtain an explicit power-law 7i-dependence 
of the energy splitting. 

Consider a mechanic system on one-dimensional circle with Hamiltonian H = Ek{p) + V{x), 
V{x + 27r) = V{x). The kinetic energy Ei^{p) satisfies E^{—p) = Ef^{p), and, for simplicity, 
we assume Ek{0) = 0, Sfc(oo) = oo and 4-Ek[p) > when p > 0. A familiar example of 



such kinetic energy is hp^- Due to the time reversal symmetry, the two classical orbits at 



dp 

H{x,p) = E > maxx V{x), and O^, one with p > and the other with p < 0, yield identical 
action integral, i.e.. 



/ pdx= (f pdx = S{E). (2.1) 
Jot Jo- 



Consequently, EBK quantization condition S{E) = 2mrh predicts a two-fold degenerate level 
E = €n- The two semi-classical eigenfunctions are given by 

*n(2;) = ^^exp[±isnix)/h], (2.2) 

where Snix) = pn{x')dx' , pn{x) > is determined by Ek{p) + V{x) = e„, Xn = ■^Ek{p)\p=p„(x) 
is the classical velocity and the normalization constant T„ = Jq^ is the period of the corre- 
sponding classical orbitf^. (The suffix "n" of e, p, x, T and s will be hereafter dropped out 
for simplicity.) 
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Of course, in general, the two levels do not exactly coincide. The difference between e and the 
exact eigenenergy is of order o{h) in the semi-classical limit (h ^ 0,n ^ oo while nh is fixed). 
In the case that V{x) is not smooth (infinitely differentiable) , we have seen in the last section 
{H2 and H4^) that the splitting of energy levels (Ae) is not exponentially small. It is therefore 
possible that a non-vanishing Ae will emerge from the higher order semi-classical corrections. If 
we are only interested in the leading term in Ae, however, variational calculation in the space 
spanned by and will give the result. We shall consider a simple case that V{x) is a C^~^ 
function and 

/\Vix)= lim ^V(x')- lim ^Vix') (2.3) 

is well-defined, which vanishes on [0, 27r] except at discrete points Xj,j = 1,...,N < 00. Then 
elementary calculations show that the energy splitting is given by (see Appendix ) 

j^k+l N Qx.u(2is(x*)/h) .k /n^ 1,11 

= S ^^^?bj7T;A.^(-l)l + o{n^-^) . AeW + oin'^-). (2.4) 

j=l P dp-^k\p=p{x*) 



Define a dimensionless measurement of ND by rjn = ^ _^i-e — Noticing that the semi-classical 
level spacing is = lirh/T and according to Eq. (2.4), we find 

fl'^ JY, exp(2/'.s-(.r*)//i) , ^ , , 

" - ^lgj4^A.V(.])l . . . (2.5, 

Example 2.1 H = + fW(x), where V^^\x) = max{cosx,0} and V^''\x) = [V^^^x)]'', 
k = 2,3,.... 

According to Eq. (2.5), when e > 1, 

2'=7r(2e)2+i n, 2 

The comparison of 77 and 77^°^ is shown in Fig.2. 

Example 2.2 H = \p\ + V^'^Xx). 
When e > 1, the semi-classical level is given by = 'nh + ak and according to Eq. (2.4) 

Ae(0) = ^^l^|sin(^ + — )| 



where 

ak = — / cos xax = . 

27ri_| 2r(i)r(| + i) 

The comparison of Ae and Ae*-'^^ is shown in Fig. 3. 
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3 Sum Over Transition Paths 



In this section we first give Eq. (2.4) a geometrical interpretation. We find the quantum 
transition between the semi-classical eigenstates can be classically described by the leaking of 
phase space points from one torus to the other via passage of non-smoothness. This picture will 
facilitate the generalization of Eq. (2.4). 

The splitting of nearly degenerate energy levels is closely related to the transition probability 
between the corresponding semi-classical eigenstates. In classical picture, describes a particle 
moving on the circle with p > 0. After one classical period, due to quantum tunneling, the 
particle has a non-zero probability to jump to the reflection wave ^~ with p < 0. Write 
exp(^)|^'+ >= c\^+ > +^|^" >. Simple calculations show that ^ < ^~\iH - e)|^'+ > 
and Ae ^|-^|- According to Eq. (A. 17), the leading term in A is the sum of contribution 
from each non-smooth point of V{x), i.e., 

N 

A w ^(0) = J2 exp(z0j), (3.1) 



with 



-,KV{x*) and 4>, = 2s{x])/n, (3.2) 



where p* = p{x*) and x* = ^-E'fc|p=p(x*)- We note that rj is exactly the reflection coefficient 
obtained by Berry 0. 

As the classical representation of and the tori and are connected by the 

straight line x = Xj where H is not smooth. We shall call the vector on x = Xj that starts 

from 0+ and ends at a transition path and denote it by (Fig. 4). Accordingly, we can say 

that — > ^~ is dominated by the tunneling along transition path(s). In fact, the reflection 

coefficient rj is determined by the local properties of 7j. Besides a constant, rj consists of 

three ingredients. /\^V{x*) can be regarded as the intensity of non-smoothness at 7j. (2p*)fc+i 

describes a power-law decay with the increase of path length 2pj = fi^^- J*, which comes 

from the product of amplitude of semi-classical wave functions, gives a classical weight of the 

transition path: the longer the particle stays in the vicinity of the non-smooth point, the more 

probably it jumps to the other torus. In contrast to rj, the phase is not determined by the 

local properties of 7^. Since only the relative phase is of physical importance, i.e., gives rise to 

^Berry's calculation was based on Ek ~ However, the result (Eq. (27) in Q) is essentially identical to Eq. 
(3.2). 
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interference effect, we find 

2 If 

where ■jjk is a closed path consists of jj, -7^ (7^ with opposite direction) and the segments of 
and 0~ (real paths) that attached at their ends (see Fig. 4). If 7jfc is contractible, (pj — (pk 
is simply the phase space area (in the unit of h) enclosed by this closed path. 

Behind the simple form of Eq. (3.2) there are two non-generic facts resulting from the 
assumption that ^E}.(j)) > when p > 0: the starting and end points of 7j are symmetric 
with respect to p = and the projection of or onto the coordinate space contains no 
singularity. Now we ignore this assumption and require only Ek{—p) = Ei.{p) to guarantee the 
time reversal symmetry. Let Aj = {x*,pj) G and A'j = (xpp'j) G 0~ be the starting and 
end points of 7^. By adopting the general semi-classical eigenfunctions corresponding to the 
tori 0+ and similar calculations as that performed in Appendix show that Eq. (3.2-3) 

should be modified as 



(p,-p;.)'^+yix(^,)iK-)l 



and 



(f>j - (pk = ^ j> pdx - MjkTT/2, (3.5) 

where is the sum of the Maslov indices of the segments of real paths on 7^-^. Having the 
contribution of each transition path, we need only to sum over all these paths to obtain the 
energy splitting Ae'-''-* or r]^^\ 

Example 3.1 H = {p^ - if + V{x), where V{x) = 1 - (f )2, |x| < vr. 
When e < 1, the Maslov index of Of (or 0~), which encircles point (vr, 1) (or (vr,— 1)), is 2 
and EBK quantization condition reads S{en) = 2(n -|- |)vr7i (see inset of Fig. 5). The straight 
line X = vr intersects Of (or 0~) at points Ai^2 (or A'^ 2) where p = (1 it £2)2 (or —(1 it €2)2). 
There exist four transition paths, i.e., 71 {Ai yl^), 72 {A2 — > ^2), 73 {Ai A'2) and 74 
{A2 — > A'l). Moreover, EBK quantization condition implies — 4>i = 0, nvr,nvr (mod 27r) for 
j = 2,3,4 respectively. According to Eq. (3.4), up to a phase, 

^(°) = A[-^ + -^ + (-l)" ' 



4vret '(l-Fet)f (l-et)t ^ ' (1 + (1 _ g)l)(i _ g)i 



When e > 1, only 71 survives so that 

^(0) 



4vre2 (1 £2 
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Numerical results show that r]^^^ = lA^^^/n is a good approximation of tj when h is sufficiently 
small and e is not too close to 1, the energy of separatrix(Fig.5). 



The same treatment can be applied to non-smooth systems where ND is originated from spa- 
tial symmetries. By substituting (p, —x) ix,p), relations (3.4-5) can be directly transformed 
to systems where the non-smoothness that results transition path depends only upon the p- 
coordinate. Specifically, consider a transition path 7-,- on the straight line p = Pj with starting 
and end points at Aj = (xj,p*) and A'^ = (x'j,p*), the corresponding reflection coefficient should 
be 

^ (i^)' /\'Hix,p*). (3.6) 

{x'^-xj)>^+^^\p{A,)p{A'^)\^^P 

The phase difference is also given by Eq. (3.5), whereas the Maslov index should count the 
singularity of the projection of torus onto the momentum space 0. Despite this similarity, 
interesting behavior may occur when the configuration space has a non-trivial topology. We 
shall demonstrate it by some examples. 

Suppose the configuration space is a circle, i.e., {x,p) and {x + 2tt,p) describe the same point. 
In this case, a path {x,p) — > {x',p) implies a family of paths {x,p) (x' + 2mT,p),n £ Z. If we 
attribute the contribution of all these paths to a representative path, say, {x,p) — > {x',p), the 
only change of Eq.(3.6) is that ^^./ _^.-,fc+i should be replaced by 

W satisfies periodic condition Wk{x, y + I) = e^'^'^^W^ix + 27r, y) = Wk{x, y). When y G [0, 1], 
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W2{x,y) = -—^[l + y{e'^-l)]e- 



■txy 



4 sin 



(3.8) 

and so on. We note that is in general not invariant under the translation {x,p) — > {x,p+6p) 
when 6p is not an integer multiple of h, which is however always a symmetric transformation in 
classical mechanics. This difference reflects the discreteness of quantum momentum space. 

Example 3.2 H = \p — pc\ + cos^ x. 
The symmetric double-well potential causes ND at e < I. According to Eq. (3.6-7), the total 
^ We use uji = pdx instead of to'i — ~xdp based on two facts. Firstly, coi = u)[ when 7jfe is contractible. 

Secondly, if the coordinate space has non-trivial topology, uj\ is well-defined while tjj is not. We find that this 
choice is justified by numerical results. 
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contribution of transition paths (on p = pc) is given by 

^(0) ^ [W2{2x,, f ) + W2i27r - 2xe, f ) + i-ir2W2{7r, f )], 

£2(1 — ej 2 li ii ii 

where Xc = cos"^ When pc = 0, A^^^ = —t^ — + When Pc = I, A^^'^ = 0. In 

fact, Ae = in this case because H is represented by the same matrix in the invariant subspaces 
spanned by functions {e*^"^}„ and {e*(^""'"^)^}„ respectively. 

Consider a spin system defined in classical and quantum mechanics by respectively { Jj, J^} = 
SjksJs and [Jj,Jk] = ihejksJs, j,k = 1,2,3. When = + Jf + Jf is fixed, the classical 
mechanics is confined within a sphere Sj. Restricting the su{2) Poisson structure to Sj yields 
a symplcctic two form uJ2 = J svnOdcj) f\ d9, where {0,(p) is the conventional sphere coordinate. 
In quantum mechanics, = j{j + j = ^, 1, An eigenspace of is associated with 

a classical sphere Sj, in which we shall assume J = {j + so that its phase space area 
(integral of uj2 on Sj ) in unit 27r/i is 2j + 1, which directly corresponds to the dimension of 
the eigenspace. In our treatment of non-smooth systems, a prerequisite is that the phase space 
is the direct product of coordinate and momentum spaces. To meet this requirement, we write 
(Jcosf? -|-po)0) = (Pi-^)) hi which lj2 = dp A dx, and regard {x,p) as the natural coordinate of 
the phase space of a mechanic system on a circle. Moreover, to ensure the right spectrum of 
J3 = p — Po, we choose po = (or ^h) in the case of j is an integer (or half integer). By this 
transformation in classical mechanics, we can treat the non-smoothness-enhanced tunneling in 
some spin systems. 



Example 3.3 

' Jl - Ji + Jl ^3>o, 
- J| J3 < 0. 



H{Jl,J2, J3) = { j2 t2 



The corresponding classical system on a circle is 



.P - {p-po)'^]cos2x + {p-po)^ p>Po, 
[J^ — {p — po)'^] cos 2x p < pq. 



Hix,p) = { ' 2 _ . __^2 



Prom phase space portrait we know that energy levels in (— J^,0) consist of 2-fold ND and 
according to Eq. (3.6-7), 

A^'^ = j2-^[^3(27r - 2x,, ^) + Ws{2x^, ^)e2^(<>-f ) + 2I^3(^, ^)e^(^-i)], 

where Xc = ^ cos~^ $, with ^ = e/J^ and = 7rJ(l — smxc)/h = TT{j -|- ^)[1 — (^^)2]. When j 
is an integer, 

1^(0) I lcos<^l 



2{j + mi-^f 
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and when j is a half integer, 

|^(0)| = 1 -I ^ + ^ 3 sin</) + -| 

(In this case, -^-y can be regarded as an effective h.) These relations give a good description of 
the energy splitting when j ^ 1 (Fig. 6). 



4 Discussion 



When the non-smooth system is controlled by a parameter A, e.g. V{x) — > XV{x), it is easy 
to obtain a zero of A^^^ when A is continuously varied. One can naturally ask that whether 
the zero of A^^^ predicts an exact degeneracy of energy level or it merely corresponds to a 
minimum of Ae. The answer turns out to be dependent upon the symmetry of the system. If 
the eigenstates involved in ND can be distinguished by different symmetries irrespective of the 
parameter, the energy difference between the two eigenstates should be a smooth function of 
A, which is approximately given by ^A^^^ or similar expression. In this case, the zero of A^^^ 
indicates a nearby exact degeneracy. Of course, because of the symmetry of H, this conclusion 
cannot be regarded as a violation of the well-known theorem of von Neumann and Wigner, which 
states that generically we must vary two parameters to create a degeneracy||8|. On the other 
hand, if the eigenstates cannot be restricted within different parameter-independent invariant 
subspaces, e.g.,when H = + cosx + A| sinx|, the zero of A^^^ generally corresponds to a 
minimum of Ae where we must take the higher order corrections into account. 
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Figure Captions 



Fig.l Splitting of nearly degenerate energy levels, (a)-(d) for H = iJi, H2, H3 and H4 respectively. 
The numerical result of Ae (open circles) , the spacing of semi-classical levels (dotted lines) and 
the semi-classical approximation of Ae (solid lines) are shown at ^ = 0.02. The insets show the 
degenerate tori (solid line) in phase-space where H is not smooth on the dotted lines. 

Fig.2 Scaled energy splitting rj (open circles) and rj^^^ (connected solid dots) in example 2.1 at fc = 1 
to 4 and n = 0.05. 

Fig.3 Energy splitting Ae (open circles) and Ae^'^^ (solid lines) in example 2.2 at A; = 1 to 4 and 
n = 0.04. 

Fig. 4 Schematic figure show transition paths 71 (^4 A'), 72 {B — > B') and closed path 721 
(B ^ B' ^ A' A' ^ B'). EBK quantization rule guarantees that (j)2 — <pi (mod 2it) is 
independent on the choice of real paths B' A' on 0~ and A —> B on Of. 

Fig. 5 rj (open circles) and rj^^^ (solid lines) in example 3.1 at ?i = 0.02. The inset shows three types 
of tori in phase space. The tori encircling point (0, 0) produce a semi-classically non-degenerate 
component of energy spectrum at 1 < e < 2, which has been excluded according to semi-classical 
criterion that the expectation value of at the corresponding eigenstates is less than unity. 

Fig.6 77 (open circles) and 77^°^ (connected dots ) in example 3.3 at (a) j = 100 and {h)j = 99^. 
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Appendix: Semi-classical Calculation of Energy Splitting 



We first consider the conventional Hamiltonian H = ^p^ + V{x). Direct calculation show that 

^"T^ + = (e + Q{x))^^{x), (Al) 

with Q = — where the prime denotes derivation with respect to x at fixed e. 
Because < ^'+1^'+ >=< >= 1 and < ^f-l^-i. >~ 0, the energy splitting calculated in 

the space spanned by \1'"'" and is given by 

Ae = 2\< > I = 2^1/0 + ^^]exp(i2.(x)/n)dx|. (A2) 

Before evaluating Ae according to Eq. (A. 2), it is helpful to recall an useful mathematical 
result on asymptotic behavior of the Fourier coefficients of a non-smooth function. Let /(x) be 
a sufficiently regular 27r-periodic function on R. How its Fourier coefficients, defined by 

f{n)= / /(x) exp(ma;)(ix, n G 2', (^-3) 

JO 

decay when n ±00 is basically determined by the analytic property of f{x). If it is smooth, 
then f{n) for large n will approach zero faster than any power of i.e., lim|„|_^o^ f{'n)\n\'^ = 

for arbitrary a > . On the other hand, if f{x) is not smooth, the decay of f{n) may follow 
a power law. In the simple case when f{x) is the union of N smooth segments on intervals 
xi < X2--- < xn+1 = xi + 27r, /(n) can be expressed by asymptotic series 

00 -i-f-i N 

hn) = i:'-nriY. eM^nx*)^J{x*). (AA) 
1=0 j=i 

Let s{x) = nh6{x), we rewrite Eq. (A. 2) as 

Noticing the integrand apart from exp(i2n^) is unchanged in the semi-classical limit, according 
to Eq. (A.4), we have 

= g p^4) A.^(^i)l + oin'-''). (A6) 

Then we consider Hamiltonian H = Ek(p) + V{x). In order to evaluate energy splitting 
according to Ae = 2| < — el^"*" > |, it is instructive to go into some details about the 
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momentum representation of Write 

'^^ = y2 (j)f\r>, < x\r >= —=exp{irx), {A.7) 

r=-oo V27r 

with 

1 f'^'^ 

(f>f = —= ^^{x)ex.p{-irx)dx. {A.8) 
V27r Jo 

Since 01^ = ^+*, we shall focus on ^'+. The semi-classical limit of Eq. (A.8) should be calculated 
in two separate cases. In the classically permissible region (CPR), where p{x) — rfi = is satisfied 
by some x G [0, 27r), stationary phase approximation can be adopted, which results 



6+ 



M exp[i(s(x„j)/^ - rxm - (TmTr/2)], {A.9) 



wherelx^} are solutions of p{x) — rh = and cr^ = signiV" (xm)) ■ When rh is beyond CPR, 
by using expansion (A.4), we find 

From Eq. (A. 9-10) wc conclude that ^'^ consists of the main part distributed within CPR and 
two powcr-law-likc long tails beyond CPR. ( As the non-smoothness of eigenfunction is resulted 
via the eigcn equation from V{x), this picture is also true for exact eigenfunction.) Furthermore, 
if the semi-classical momentum representation of is calculated in the similar procedure, 

one can find that the main part of within CPR but its long tails approximately satisfies 
eigen equation {E^ip) + V{x))'^ = e^, i.e., 

OO 

m=—oo 

when rh G CPR, where 

i'^+i ^ k 
Vm =< 0\V\m >« ^-^Y,exp{imx*)/\V{x*) (|m| oo). {A.12) 

Based on the above discussion, we know that 

OO 

<^-\Ekip)\^+ > = ct^trMrmt- E ^trEk{rh)cj>+ 

r=-oo |r?i|GCPR 



( E E + E E )^-uies.,, - v^..)^; 

i/?ieCPR/i=-oo -ij.neCPRv=-oo 



{A.13) 
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Compare the last expression with 

oo 

The main contribution of Eq. (A. 14) consists of three parts come from regions, {l)fih, uU G CPR, 
{2)—fih, —ufi G CPR and (3)/Lt^, —uU G CPR respectively. Eq. (A. 13) contains only the former 
two parts while we can screen the last contribution by making a high frequency cut off of V(x), 
i.e., replacing it by 

y(o)(ar)= ^ y^exp(-zmx), (A15) 

\m\<kc 

where kc is a large but fixed integer so that F(i)(x) = V{x) - is negligibly small. 

Therefore, < ^'~|Efe(p)|^'+ >Ri< *~|e - >, and consequently 

1 /■^'^ V^^^ 

< *-|ii"-e|*+ >f5i< >= - / ^ e-x.^{2is{x)/h)dx. (A16) 

Observing that « and /\iV^^\x) = /\iV{x) for arbitrary x G [0,27r) and j > 0, by 

partial integrating Eq. (A. 16) for successive k + 1 times we obtain 

which immediately leads to Eq. (2.4). 

Finally, it is worth pointing out that although the exact eigenstates have power-law tails beyond 
CPR, the leading term of Ae actually does not relies on this detail. In fact, Eq. (A. 16) essentially 
equals to 

which is in nature controlled by the power-law decay of {V^} but {(f)^}- Therefore, Eq. 
(A. 16) can be reproduced from the highly localized semi-classical eigenfunctions correspond- 
ing to smoothed Hamiltonian H^^^ = E^{jp) + V^'^\x). 
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